Estimate of the Fundamental Solution for Parabolic 



Operators with Discontinuous Coefficients 
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Abstract. We will show that the same type of estimates known for the fundamental 
solutions for scalar parabolic equations with smooth enough coefficients hold for the 
first order derivatives of fundamental solution with respect to space variables of scalar 
parabolic equations of divergence form with discontinuous coefficients. The estimate 
is very important for many applications. For example, it is important for the inverse 
problem identifying inclusions inside a heat conductive medium from boundary mea- 
surements. 

Mathematics Subject Classification(2000): 35R30. 

1 Introduction 

Let X be a parabolic operator of the form 

£ = dt-V-AV (1.1) 

with an n X n matrix A = (aij) e L'^{D) and a bounded domain D c R" with boundary dD of 
Lipschitz class. D has compactly embedded subdomain D,„ (m = 1, 2, • ■ • , L) with boundaries 
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dD„ (m = 1, 2, • • • ,L) of C^'" class such that D( n D„, = %{i m), D„^ a D {\ < m < L), where 
< or < 1. Moreover, there exists a constant 6 > such that 

n n 

J] aijixMj > 5 2 ^2 e r , a.e. x e D). (1.2) 

Also, we assume that A e C^iDm) (1 < m < L) with < // < 1. 

We want to show the following theorem which is quite important in inverse problems for heat 
equations with discontinuous coefficients. (See the remark after the following theorem.) 

Theorem 1.1. There exists a fundamental solution Y{x, t;y, s) ofX. with the estimates 

0<Y{x,t-y,s)<-^^-^e-'-^^, (1.3) 

^ (,_,)(n.l)/2 ^"^ (1-4) 

for any t, s e'K, t > s and almost every x,y e D, where C > is a constant depending only on A 
and n. 

Remark 1.2. 

( i) For the simplicity of notations, we confined our argument to scalar parabolic operators of 
divergence form without zeroth order term. However, our argument can be generalized not only to 
more general scalar parabolic operators, but also to parabolic systems. 

( ii) The estimate (|1.3I) is the well known estimate of the fundamental solution ([2)]). We will call the 
estimate (11.41) gradient estimate. This gradient estimate is quite crucial for the dynamical probe 
method ^ and stability estimate for the inverse boundary value problem /12]/- Here, the dynamical 
probe method is a reconstruction scheme for the inverse boundary value problem identifying an 
unknown discontinuities of a medium inside a known heat conductor from our measurements on 
the boundary which are so called Dirichlet to Neumann map or the Neumann to Dirichlet map. The 
graphs of these maps are nothing but the set of infinitely many Cauchy data of the solutions to the 
forward problem of this inverse problem. Also, the stability estimate is the estimate of continuity 
of the map which maps the set of Cauchy data to the unknown this continuity of media. 

We will show later that the gradient estimate of the fundamental solution follows from the 
following interior estimate following the argument given in (31. 
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Theorem 1.3 (Main Theorem). Let < r < T, rQ c D and u e W{rQ x (-r^, r^)); u e 
L^{{-r'^,r^);Hl^{ra)); d^u 6 L'^ii-r^,r^y, H-\rQ.)) be a solution of {d, - V ■ AV)u = in rQ. x 
(-r^, r^). Then, there exists a constant c > such that for any < p < r/2 anJ {x, t) e 
(r - 2p)Q. X (-r^ + 4p^, r^), we have 

c 

xU\\L'^{pO.{x)x(-p^^t,t)) ^ ^^'^^^L^(2pn(x)x(-4p^+t,t))^ (1-5) 

where Q.{x) := {y = (ji, ■ ■ ■ e R"; < 1 (1 < ? < n)} with x = (xi, • • • ,Xn) andO. := f2(0). 

Remark 1.4. In ?/zere a/50 a statement of the corresponding main theorem. Here we will 
give a proof of the main theorem based on the argument of ^ for second order elliptic systems of 
divergence forms which could provide a proof of the main theorem even in the case the inclusions 
touch. 

By translation, rotation and scaling, it is enough to prove the following. 

Proposition 1.5. Let A as in Theorem \L3\ Then, for any solution u e W{Q. x (-1, 1)) o/(t?, - V ■ 
AV)u = m Q X (-1, 1), there exists a constant c > depending only on the bound of A and n such 
that 

|V,m(0,t)| <c||M||^2(inx(|r4,r)) (t £ 1)) (1-6) 

whenever VvM(0, t) exists. 

The rest of this paper is organized as follows. In the next section we will give the proof of 
Proposition [T3]by assuming the existence of a Green function in a two layered cube with Dirichlet 
boundary condition and its estimates. The proof of the gradient estimate of a fundamental solution 
is given in section 3. Finally, in Appendix, we give a construction and estimates of the Green 
function in the two layered cube with Dirichlet boundary condition. 



2 Proof of Proposition [13 

We will adapt the proof of Li-Nireberg's paper [|7]|. To begin with, we note that Lemma 4.3 in |I71 
also holds if || • ||yi+£,',2 is replaced by || • Hyi+a'.,. for any p > \. That is 
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Lemma 2.1 (Lemma 4.3'). For < a' < min{/i, ^^^^^^ }, r/zere exists a constant E > depending 
only on ||A||c<i' /or these D^'s which intersect with Q such that 

\\A-A\\yu.',, <E, (2.1) 

A is defined as in [T] for the very special case i.e. the two layered cube D. 

Proof. This can be easily proved by checking the proof in [7J. □ 

Remark 2.2. Hereafter, names of theorems in parentheses correspond to those of the theorem in 
/[7|/. For example, Lemma \2J\ Lemma 4.3') correspond to Lemma 4.3 in /[7|/. 

Next we generalize Lemma 3.1 in tTJ in a special way to the parabolic operator (|l.ll) . That is 
we have the following: 

Lemma 2.3 (Lemma 3.1'). For < s < I, let the previous A satisfy 

||A < e, (2.2) 

where p will be specified in the proof. Here, we note that we have properly scaled A so that by 
applying Lemma IZTI ( Lemma 4.3'), (12.21 ) is satisfied. Then, for any f e L^((-l, r); L^(Q)) and 
solution u e W{Q. x (-1, r)) of 

{dt-V ■AV)u = V ■ f in fix(-l,T), 

there exists a solution v e V7(|Q x (^r - ^, r)) of 

{d,-V-AV)v = in Iqx(^t-^,t) 

4 Id Id 

and I < (To < 1 such that for any < cr < ctq 

\\U - v||L2(o-nx((l-o-2)T-o-2,r)) ^ C'(||/||t2(nx(-l,r)) + £^^^\\u\\lHQ.x{-1,t))), (2.3) 

where C > depends only on n and A. 

Proof. It is enough to prove the estimate for the case that Q. is divided into two parts by the 
boundary dD,„ of one of D,„ (1 < m < L) such that the center e dD^ of CI and Q does not contain 
any portion of dD{ {£ m, 1 < £ < m). 



Take a cutoff function ^ G Cq(Q. x (-1, 2t)). By the definition of weak solution, we have 

T r 

J [id,u, if) + (AVm, V(^)) J? = - V(^)rf? 



for any e H^^iO.). 

■2. 



If we take (p = ^ u, then we have 

T r 

/((.,„,^^„,.(AV„,V(^^„,,).,= -/(/,V(,^„)W,. 

-1 -1 

Here, we have 

LHS>i J (fV)(-,T)djc + 

€1 -in 

T T 

5 ^'^\^u\^dxdt -C J J (^\Vu\)(u\V^\)dxdt, 



+ 

-1 a -in 
where C > denotes any general constant. We further have 



r T 

^ u^dxdT, 



-in -in 
where C(6) > denotes any general constant depending on 6. 
On the other hand, we have 

T 

RHS =-//(,/. «V„).„/.VfV^. 

-1 n 

T T 

-in -in 

Hence, 

T T 

-in -in 
By |V(^m)|2 < 2^2|v„|2 + 2u^\V^\^, we have 

r r 



-in -in 



Then, combining this with Poincare's inequality, we have 



< C(5)(ll"llL2(nx(-l,r)) + ll/llL2(nx(-l,r))). 



If we let ^ to satisfy ^ = 1 in a neighborhood of |r2 x (^t - ^, r), we have 

I|m|Il2((|^_16 .,).^i(4^)) < C((5)(||M||i2(nx(-l,r)) + ll/llL2(nx(-l,T)))- 

By (12.41 ) and the Fubini theorem, there exists | < ctq < 1 such that 

ll"llL2(((l-o-2)T-o-2,r);//l('5(o-on))) ^ C'(5)(ll'^llL2(nx(-l,r)) + ll/llL2(nx(-l,r))) 

Let V e Wio-Q. x (-1, r)) be the solution to 

'(5, -V-AV)v = in o-qQ x ((1 - a-^)T - cr^, r). 



V = u 



V = u 



on 5(0-0^^) X ((1 - crl)T - ctq, r), . 
on o-o0.x [(1 - crl)T - crh. 



Then, w := a - v satisfies 

(5,-V-IV)w = V-((A-I)Vm) + V-/ in o-qQ x ((1 - (r^)r - trg, r). 



w = 
w = 



on dio-oQ.) X ((1 - crl)T - (Tq, t), 
on o-qQ X {(1 - crQ)T - ctq}. 



(2.4) 



Further, let G*(x, t; y, s) be the Green function such that for y e crQ and s eR. That is, G*(x, t; y, s) 
is a distribution which satisfies 

(d, + V- AV)G*(x, t; y,s) = in o-qO. x (-oo, 
G*(x, t; y,s) = on (9(croQ) x (-oo, s), 

lim r G*ix,t-y,sMx)dx = <p{y) (if e C^icron)). 
J 

(ToQ. 

Later in the appendix, we will provide the construction of G*(x, t\y, s) and proof of its estimate: 

\d"G*{x,t;y,s)\<Cais-t)~'^e~^^ (?, 5 e R, ? < 5, a.e. x, j e o-qQ). (2.5) 



for any a G Z" , |a| < 1 with some constant Cq, > 0. 



By the Green formula, we have 



w{x,t) = - J J [({A-A)Vu)(y,s) + f(y,s)}-VyG*{y,s;x,t)dyds. (2.6) 

for X 6 o-qO., t e ((I - crl)T - a^, r). 

To proceed further, we need the following Lemma [Z4l which is well-known in Fourier analysis. 

Lemma 2.4. Let {Xi,Mi,mi) {i = 1,2) be cr finite complete measure spaces and {X\ x X2,M\ 
M2, mi X mi) be the product space with complete measure m\ x m2. Also, let p\, p2 and q e [1, 00) 
satisfy \ /p2 + I /q = 1 / pi + 1 and a measurable function K{xi , X2) on X[ x X2 satisfy 



/ 
I 

X2 



\K(xi,X2Wmi{dxi) < Li (a.e. X2 e X2), 



\K(xi,X2)\''m2{dx2) < L2 {a.e. xi e X\). 



Then, for any f(x2) e Lp^{X2), we have 



\\Kf\\ 



LP\(X,)<l^i 1^2 



LP2{X2), 



where 



{Kf)ixi) = J K(xuX2)fix2)m2(dx2). 

We apply this Lemma[2]4]to the kernel (A - A){y, s)VyG*{y, s; x, t) {{x, t), {y, s) e o-qO. x ((1 - 
(tI)t - cFq, t), s < t)hy taking q = I, p^ = p2 = 2 and Xi = X2 = o-qQ. x ((1 - crl)T - cr^, r). 
Let 

hiy,s)= J J \{A-A)(y,s)VyG*iy,s;x,t)\dxdt, 

{\-o-l)T-o-l o-on 

and 

T 

h(x, t)= J J ^^^-^^ s)VyG*{y, s; x, t) dyds. 
Since A, A are bounded and Jjc = 0(St - 5)"^^) {s < t), 

T 

h(y, s)<C Jit- sT^i^dt = 2C(t - sf'^ < C. 



Similarly, by taking p > n + 2, 



h(x, t) 



A i I 

(l-o-2)r-o-2 o-o" 



(A - A)(y, s) dyds 



n / / 



VyG (y, s; x, t) dyds 



Up' 



< C(T - (1 - crl)T + criy'P\\A - A||„(a)( j (t - ^)"/2-("+i)P*/2j,j 



with l/p + lip* = 1. Here, by p > n + 2, n/2 - (n + l)p*/2 > -1. Hence, hixj) < C||A - A||iP(n)- 
Therefore, by Lemma [Z4l and (12.41) . 



|l^llL2((ronx((l-tr2)r-^2,r)) ^ " ^llll?n)l|VM|lL2(cronx((l-^2)r-(r2,r)) + 



L2(cronx((l-cr; 



■g)r-.r2,r)))- 



Since ||A - A||yi+„'.p < £ implies ||A - A||LP(n) < Ce, we have 



lkllL2(o-onx((l-o-2)r-o-2,T)) ^ C(£^''^||M||i2(nx(-l,T)) + ll/llL2(nx(-l,r)))- 

Further, by applying (|2.4I) to w, we have for a smaller ctq 

\\w\\LH({l-o-^)r-o-^,ry,HH<Ta)) ^ C(£^^^||a||i2(nx(-l,r)) + ll/llL2(nx(-l,r))) (0 < (T < (Tq). 

This ends the proof. 



□ 



Proof of (11.61) . We adapt the proof of Proposition 4.1 in Q to our case. 

Let o-k = ^,o=k = ^,o-k= ^FTT = 0, 1,2, • • •) and M = l|M|lL2(inx(|r-i,r))- We will prove 
by induction that there exist wt e W(o-kO. x ((1 - cr^)T - cr^, r)) (fc = 0, 1, 2, ■ ■ ■ ) which satisfy 



(dt-V ■AV)wk = in (Tknx{(l-crhT-o-l,T), 



k(n+4+2a') 



\\Wk\\L2iTf,nxiil-?rl)r-?rlr)) ^ CMA' ^ , 
WxWk\\L'-ia;£ix((l~d-l)T^alT)) ^ CMA'''" , 



(4-3)1 

(4.4);,, 

(4-4)y 
(4-5); 

L^{at0.y.{il-^)r-al,r)) 

Before starting the induction argument, we note that for any sq > 0, we have ||A - A||yi+a'.p < sq 
by considering a dilated Q. instead of Q., and ||A - A||^p(n) < ||A - AHyi+c,',,,. 



and 



< MA 5 . 
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Since u e W(\Q. x (|t - |,t)) solves (5, - V • AV)u = in x (|t - |,t), we have from 
Lemma [X?] that there exists a solution wq e W(o-oQ. x ((1 - cr^)T - cr^,T)) of (5, - V ■ AV)w^: = 
in ctqQ, X ((1 - (tI)t - cr^, r) with the estimates: 

II" - >volL2(^„^,,((i„^2)^_-.^^)) < Coe^/'M < M4-^ (4.5)', 

by taking so > small enough to satisfy Cqs^^^ < 1, e^"* < 4'"*'^^^''' and hence 

INolL2(^„nx((i-^^)r-^^)) ^ 2M, (4.4); 

and from the interior estimate for 5^ - V • AV, we have 

l|V;t>Vo|lL2(5:nx((l-5^2)r-5^2 r)) < CqM. 

Hereafter, Co > is a general constant for the estimate of solutions for our parabolic operators 
which is independent of the general constant C in the estimates (4.4)'-(4.5)'. 

Suppose (4.3)'-(4.5)' hold up to ^ > 0. Then, we will prove them for k + I. Let 

k 

W(x, t) = (u-J^ Wj^^kX, alt + (1 - d-l)T) {{x, t)eax (-1, r)), 

;=o 

A^^i{x) =A(d-kX), A^+i(x) =A(d-kX) 

and 

k 

fk+iix, t) = aoiAk+i - Ai+i)(jc) ^ VwjiakX, alt + (1 - al)T). 

7=0 

Then, it is not hard to see that (5( - V • A^+i V)VK = V • f^+i in x (-1, r). Further, we have 

k k 

Yj \(y^Wj)(5-kX, alt + (1 - al)r)\ < CmY, < ((-^^ t)eQx{-l,T)) 

;=0 7=0 

by (4.4)^ „ and 

l|W|L2(i,,(_i,,)) < M4-^'^'^^'-'^'^ 

by (4.5)^. Observe that 

||A,,i - < 4-'^'^'^"'\\A-A\\yu.,2 < 4-(^+i)'^'||A - < 4-^'-''^''' so- 

Together with this and (4.4)2 ^^^^ 

IIA.illL2(nx(-i,.)) < CM4-('^->Xi-')^oi;4->' < j^so. 

7=0 



By Lemma |2.3l there exists a solution vm e W((ToQ.x{{l-crl)T-crl,T)) of (df-V ■Ai,+iV)vk+i =0 
in o-qQ. X ((1 - cr*')r - cr^, r) with the estimate 



C 



,1/2 



Let Wk+i(x, t) = vt+ii^^^x, aft + (1 - c^f)T) ((jc, t) e crt+i^ x ((1 - cr^+i)T - cr^+i, t)). Then, 
it is easy to see that (d, - V • AV)wk+i = in (Tk+i^ x ((1 - cr^+i)T - cr^^j,T)). Further, by 
Vk+i(x, t) = Wk+x(akX, alt + (1 - ?^)t) {{x, t) e o-qQ x ((1 - cr^)T - cr^, r)), 

yt+l 

(W^ - Vi+l)(x, = (m - ^ C?^? + (1 - C?^)t), 

7=0 



I.e. 



k+i 



(u'Yj ^;)(^' t) = (W- Vk^iXo-fx, crft + (1 - erf )t). 

Since ?o^^x((l-a'2)^_5:2^^^ ^ a=orix((l-a=o)T-a=o, r) and (x, e o=^+iQx((1-o=2^j)t-?2^i,t) 

which equivalents to (cr^'jc, crft + (1 - crf)T) & ctq^ x ((1 - crl)T - cTq, t), we have 

/t+i 

_ (n+2)/2 



j=0 



L2(5^i.+inx((l-5^2^,)r-5^2^j,r)) 



and hence by the estimate of \\W - Vk+iWLUj^^ii^^^i.-^i^r-Trlr))' we have 



k+i 



c 



i=0 



C 



(t+2)(n+4+2g') 
2 



<2CoM4/^max(^-^,l)4- 

Therefore, by taking sq > small enough to satisfy ICqsI'^ max ( , ij < 1, we have (4.5)^^j. 
By the estimate of - V/t+i||^2(^gQx({i-o=5)T-o=^ r)) ^^^^ '^h^ interior estimate for 5, - V • AV, 

1 1 VV/:+ 1 1 |ioo(-|,fjx((l-5=2)r-5:2,r)) 
< Co||V^+l ||^2(5=(,£|><((i_g=2j^_-2_^^j 

C()(l|M^llL2(Jfonx((l-Jfg)r-Jfg,r)) + ' ^^+1 llL2(^onx((l-Jfg)r-Zfg,r))) 



< 



< CoM 



C 



V So + e. 



l+2Co6:yVax(^-^,l) 
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Hence, if we adjust C and sq to satisfy C > Co[l + ICqSq^ max( ^_^_„, , 1)], then we have (4.4)2^_^i. 
Finally, to see (4.4)^ we have from the above estimate of l|v^;+i||^2(5:(,j^x((i-5^2)^_^2 and (x, t) G 
au+iO. X ((1 - o^^+i)t - a\j^i,T) which is equivalent to (a^^x,o-l^t + (1 - e o^qO. x ((1 - 

a^o)T - o^Q, r). we have 



<CqMA 



-(yt+lXl+a')~(n+2)/2 



l+2Co^;^^max(^-^,l) 



(t+l)(n+4+2a') 



l+2Co4/Vax(^-^,l) 



Therefore, if we further adjust C and sq to satisfy C > Cq{\ + 2CqEq^ max( , 1)], then we have 
(4-4)'i,fc+i. Thus, we have proven (4.3)' - (4.5)'. 

Let C be a general constant which is different from the general constant C in (4.4)' - (4.5)'. As 
an easy consequence of (4.4)2 k' ^^^^ 



(4.6)' 



Together with this and (4.4)2, 

k oo k ~ 

I ^ Wj(x, t)-J] r)\ <CM 4~^"'|(x, f - t)| + CM J] 4" 

<CM\(x, t-T)\ + CM4-^^i+"'\ 



7=0 



■KM) 



=k+l 



Hence, we have 



CO 



L2(?tnx((l-?|)T-5^,r)) 



< 



k k oo 

M - > W; + / Wj—} >V,(0,t) 

7=0 * * 7=0 7=0 

T 



L2(?4nx((l-5^2)T-5^2T)) 



(t+l)(n+4+2ii') 

< M4 2 + 



Cm| J J {\x\^ + (t - rf + 4-^''^^*'''^)dxdt 



1/2 



Hence, by {t - < a^iT + 1)^ = A'^^^^^Xt + 1), we can absorb {t - rf into Since 

C UPJX < C4~^*"^^^^''"^^^ C 4~^*^^"^'^'^Jx < C4-2*(l+a')-(*+l)n ^nd 4-2*(l+«')-(^+l)« < 4-(fc+l)(n+2) fQj. 

large ^, we have 



£vyXO,r) 



7=0 



L2(?iax((l-?2)r-52,T)) 



(t+l)(n+2) 

< CMA —au. 



Therefore, m(0, t) = ^7(0, r) and |V;,m(0, t)| < CM if V;,m(0, t) exists. 
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3 Gradient Estimate of Fundamental Solution 



In this section, as we already mentioned in the introduction, we will give an estimate of VJ^{x, t; y, s) 
for a fundamental solution Y{x, t; y, s) of the operator X as an application of our main theorem 
(Theorem 1 1.31) by following the argument given in For the readers' convenience, we repeat the 
argument. 

It is well known that there exists a fundamental solution Y(x, t; y, s) with the estimate 



Y(x, t;y, s) < jj-j-^ — -"[^^ ^Xu,c<,] (t,seR,t> s, a.e. x,y e D), (3.1) 



which is positive for t > s, where C > is a constant which depends only on A, n and X[s.oo) is the 
characteristic function of [s, oo). (See dlj.) 
Now we state the estimate of VxY{x, t; y, s). 

Proposition 3.1. Let Y{x, t;y, s) be the previous fundamental solution of the operator 5^ - V • AV. 
There exists a constant C > depending only on A and n such that 

|V,r(x, t; y, s)\ < —e-^ , (3.2) 

(t - s)^ 

for any t,s eR, t > s and almost every x,y G D. 

Remark 3.2. We recall that a fundamental solution G*{x, t, y, s) of the operator dt + V- AV can be 

given by 

G*(x,t;y,s) = Y(y,s;x,t) {{x,t),(y, s) e Q := D xR, (x,t) (y, s)). (3.3) 
Hence, estimates similar to (|3.1I) and (|3.2I) hold for G*(x, t;y, s). 

Before proving Proposition l3.1l we give the following estimate which is necessary for the proof. 

Proposition 3.3. Let QpixQ,to) = 5p(jco) x (to - p^Jo), Bp(xo) := {x e R" ; \x - xq\ < p}. There 
exists a constant C > depending only on A and n such that the following inequality holds. 

\Y(x, t; ^, rfdxdt < C—^--—e (r < t^, a.e. ^ e D), (3.4) 

J (^0 - t)"-i 

where p = |[Uo - + to - t]^'^. 
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Proof. From the inequality (13.11 ) we have 

J \Y{x,t-^,Tfdxdt<Cy J -^j^-^^e'^^X[r,^co)dxdt, (3.5) 

where Ci > is a constant depending only on A and n. In what follows we denote by / the integral 
at the right-hand side of (13.51) . We distinguish two cases 

i) to-p^ <T< to, 

ii) T <to- p^. 

Let us consider case i). It is easy to see that there exists a constant C > such that 

C-'p<\x-^\<Cp (xeB,(xo)). (3.6) 

By (|3.6I) we have 

to-T 

(3.7, 



where c„ > is a constant depending only on n and C2 > is a constant depending only on A and 

2 _ P~ 2 

n. Now we assume < to - t < Since s'"e "^2' is an increasing function in (0, ^), we have 

CnP" - 

I < — — -e (3.8) 

(?o-T)«-i 



by dXV 



Further, if we assume < ?o - r < p^, then 



nC2 

1 
P- 

c 



p~"I J s~"e ^^ds < 



{to-ry-' 





g C2(fo- 



due to the equivalence of ? - r and p^, where c„ and C2 are the same kind of constants as before. 
Hence, by the last inequality and (13.81) . we have the Proposition in case i). 
Let us consider case ii). It is easy to see that 



6p^ <\x-^\^ + t-T < 60p^ (3.9) 



for every {x, t) e Qp{xo, to). Moreover, denoting 



g 2Ci(r-T) 

Mp = max \ Jf—^'' ^ 2p(^o, k) \ 
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and taking into account (13.91 ). we have 



(3.10) 



where C > is a constant depending only on n. Now, since r < to- p^, we have 



to-T 



< 16. 



(3.11) 



Therefore, by (13.101 ) and (13.1 II ). we have the Proposition in case ii) as well. 



□ 



Proof of Propositio ns. 1\ By applying our main theorem to the function r(-, •; ^, t), we have 

ll/2 



c 



||Vr(-, ■;^,T)||L<»(Q^(;,o,fo)) < 

p 2 



\Y{x, t; ^,t)\ dxdt 



(3.12) 



22p(j;o,fo) 

Further, applying Proposition |3.3| to the right-hand side of (13.121 ) we have 



C 



\WT(-,-',^,T)\\L--'(Qp{xo,to)) ^ ~7t4 

Then, we immediately have (13.21) . because 



(?0-T)«-l 



g ato-T) 



1/2 



1 C 
— < . 

p2 - T 



□ 



4 Appendix A : Construction of Green Function in two Layered 
Cube 

In this section we will construct the Green function G*(x, t; y, s) of our operator 5, + V • AV in 
ctqQ X R with Dirichlet boundary condition on d((TQQ.) x R. If G{x, t; y, s) is the Green function of 
the operator £. = d,-V -AV in (TqQ. x R with Dirichlet boundary condition on (^(c^o^^) x R, we have 



G*(x, t;y, s) = G(y,s;x, t). 



(4.1) 



Hence, it is enough to construct the Green function G{x, t; y, s). 
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First we construct a fundamental solution T(x, t; y, s) of X. We divide the construction into two 
cases. They are y„ > and y„ < 0. We first consider the case y„ > 0. Let A, B be positive definite 

symmetric constant matrices. A = {aij)\<ij<n,B = (bij)i<ij<„. Define A = A + (B - A)x-i^), where 



0, ^„ > 0, 

1, ^n<0. 

Let r(jc, t; y, s) be the fundamental solution for dt-V • (AV;^), that is, 

dtT(x, t;y,s)-V- (AW(x, t; y, s)) = Six - y)S(t - s). 



(4.2) 



Note that Y(x, t; y, s) is also the fundamental of the Cauchy problem at t = s for the operator 
5, - V • AV. Let f be the Laplace transform of F with respect to t, that is. 







(4.3) 



Then f satisfies 



tF(x, - V • (AVJTix, 0) = Six - y)e- 
Now, we denote F for different regions as follows: 



(4.4) 



F = 



F^i for x„ > y„, 

for y„>x„> 0, 
F^ for > x„. 



(4.5) 



For cp G C^(R"), we have 
= J [Tt(p + A Vf • V(p - e^'S(x - y)(p]dx 

R" 

= J [TTif - e-^'S(x - y)(p]dx 

R" 

+ J AVf" ■Vcp+ J AVfi^ • + J BVt^ ■ V(p 

X„>yn yn>Xn>Q 0>x„ 

= J TYcpdx- J e-^'S(x' -y')(pdx' - J AVt^^ ■ e„(pdx' - J V • (AVf")^djc 



Xn=yn 



Xn=yn 



Xn>yn 
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+ 



J AVf 12 • e„<pdx' - J AVf 12 • e^cpdx' -J V • (AVf^^)^ 

n —yn 

J 5Vf2 • e„(pdx' - ^ • (BW^)(p 



X„=0 yn>Xn>0 



Xn=yn 

+ J BW^-e„ 

x„ =0 0>x„ 

= J rtcpdx -J V(A • Vfii)v? - J V • (AVfi2)^ ~ / ^ ' ^^^^^^^ 



>'„>j:„>0 



0>Xn 



+ 



Xfi —yn Xfi —0 

Therefore, we have the following transmission problem 

V • (AVfii) - rf" =0 in Xn>y„ 

V • (AVfi2) - tP^ = in );„ > jc„ > 

V ■{BVP)-tP = Q in 0>x„ 
fii-fi2 = on Xn = yn 
AV(fii-fi2).e„ = -e-"c5(x'-/) on x„ = };„ 
fi2 - f2 = on x„ = 
AVfi2 . en - BVP • e„ = on x„ = 0, 

where e„ = (0, • • • , 0, 1). Let ^^i^^.z ^j^g Fourier transforms of f^^'i^-^ for x' - (xi, • • • , x„_i). 
From now on, we use ^' = (^i, • • • ,^n-\) to denote the Fourier variable associated with x'. Then, 
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we have 



dxl ^ 4]^ 


-(if ■^'+T)0l' =0 


in {x„ 


> yn]. 






in {y„ > Xn 


>0}, 






in {xn 


<0}, 






on {x„ 
on {x„ 


= yn]. 






on [Xn 


= 0}, 




= on {x„ 


= 0}, 



(4.6) 



where A = (aij)i<ij<n-i, B = (bij)i<ij<n-i. In addition, we put another conditions 

lim 0" = 0, lim 0^ = q. (4.7) 
For simplicity of notations, let us put 

@A = Van,m' ■ r + T) - a2]i/2, 0^ = [Z,„„(5^' . f + T) - 
where the real parts of ©a and are positive. From the first three difl'erential equations in (14.61) . 



we have 



= Ci exp [ 
= C2exp[- 
0^ = C4 exp 



-za - ©A 



-/a - ©A 



-/I j J 

x„J + C3 exp [ Xn\, 



-ib + Qb 



Conditions on x„ = j„ and jc„ = imply that 
Ci-C2-C3exp[^y„]=0, 

(Ci - C2)(-/a - ©a) - C3(-/a + ©a) exp If^yJ = -e-"->' r exp [^jj, 
C2 + C3 - C4 = 0, 

C2©A - C3©a + C4©B = 0. 
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Then, we have 
Hence, we have 

&A - &B -Ts-if-r \ -ia-&A id - ®A 1 

+ 'y f exp [ x„ + y^], 

,2 1 -Ts-iV-e \-ib + @B ia-@A ^ 
©A + 0B o„„ a„„ J 

Therefore, we have the following forms for Y 

(T+ioo 

r2 cr-(oo 
cr+/oo 

r2 cr-ioo 

r2(x,?;y,5) = -i- r /^'-^'>^' f e'^'-'^V\xn,yn,^',T)dTd^', 
(2n)"i J J 

R2 cr-ioo 

where cr > and 

V (x„,yn, ^ ,t) = — - exp [ x„ + y„J 

0A - 0B r -ia - 0A j« - 0A 1 

20a(0a + 0b) L 
,,12, ^/ , 0A-0B r-/a-0A /a-0A 1 

20a(0A + 0b) ann ^ 

1 r -ia + 0A ia - 0a i 
— exp [ Xn + yn\. 



+ 

20A O-nn O-nn 
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V\x„,y„,^',T) = 



1 r-ib + &B ia-&A 

— — exp ; X„ H 

0A + 0B L b„n ttnn 



'4 



Next we consider the case jn < 0. We denote T for different regions as follows. That is 

for Xn > 0, 

r = <r2i for 0>jc„>3;„, (4.8) 
for y„ > x„. 

Let f ^'^^'^^ be the Laplace transform of r^'^^'^^ with respect to t and (p^'^^'^^ be the Fourier transforms 
of f '^''^^ with respect to x' = (xi, • • • ,jc„_i). Here we used the notation r^'^^'^^ for example to 
represent one of r\ F^^. 

Then, by a similar argument as we did for the case y„ > 0, we have 



0i = — i— e-"-'>'f'exp[^ 



-ia - ©A ib + ©B 



2© 



+ 



©«-© 



■'nn 



^22 



e " ^ exp [ ; Xn + 

©A 

+ ©a) 
u ^_ 

20b(0b + 



nn 



1 



20b(©b + ©a) 
©B -0A 



''■^'exp[: 



x„ + 



-iZ? + ©B + ©B 

A -rs-iv-e \-ib + &B , ib + &B 1 
— -e f exp — Xn + — yn\ 

©a) ^ bnn bnn ^ 



+ 



1 



2©. 



-Ts-iy 



exp 



-/Z? + ©B ib - ©B 



■Jn 



Therefore, we have the following forms for F 



(T+IOO 



(T-ica 

(T+ioo 



cr-(oo 

(T+!00 



F22(;c,r;3;, 5) 

where cr > and 



R2 cr-ioo 



,,1, , 1 r-ja-0A ib + @B 1 
y (Xn,yn, ^ ,t)= exp ;c„ + — yn \, 

0A + 0B ann bn„ ^ 

^ (^n,yn,^,T) = — — exp — x„ + — y„\ 

2©R L ^ ^ J 
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©B - ©A r -ib + ©B ib + ©s 

+ 



exp [— Xn + — yn\. 



2©b(©b + ©a) 

©B - ©A r -ib + ©B «^ + ©B 



1 r-/Z7 + ©B ib-@B ^ 



2©B Z7„„ Z7„„ 

Next we will show how to construct the Green function G{x, t; y, s) from Y{x, t; y, s) by using 
the argument given in [jEl. For example, consider a face xi = -ctq of ctqQ,. For the simplicity of 
notations, we introduce f(xi,x",t;y,s) = r(xi - cro,x",t;yi - cro,y",s) with x" = {x2,-- - ,x„)- 
Then, f solves 

{dt-V ■ AV)f = in R" x (s, oo) 
- (4.9) 
lim^, X,„ f t; y, s)(t>{y)dy = 0(^i - ctq, x") (0 e ^"'(R")) 

Let us distinguish A here by denoting it by A. Now, we extend A = {a^^) 'mx\ > denoted by A+ to 
< as follows. That is we define A_ = (a^r.) by = a|j, ar. = at (2 < /, j < n), = -a[j (2 < 
J < n). Then, if we define t'{x,t;y,s) by f'Cxi, jc", f; 5) = r(+xi, jc", f;};, 5) > 0), then F 
satisfies 

{dt-V ■ AV)f ' = in R" x {s, 00) 

(4.10) 

lim^, r(x, J, ^)<^(j)^i?y = 0(^1 - cTo, X-) ((f> 6 C~({xi > 0} x R""') 
and f '(x, t; -y, s) = t'(-xi , x", t; y, s) (t > s, a.e. x, j e R"). Hence, T'(x, t; y, s) - Y'{x, -ji , / ', s) 
for x\,y\ > is the Green function in the domain {^i > 0} satisfying the Dirichlet boundary 
condition on xi = 0. Repeating this argument for other faces of ctqQ., we can construct the Green 
function G(x, t; y, s) for x,ye a^Q.. It is clear from its construction that G{x, t; y, s) satisfies the 
estimate (|2.5I) and by (|4.1I) . G*(x, t; y, s) also satisfies the same estimate. 



5 Appendix B: Estimate of the Green Function 

In order to give its meaning to the fundamental solution r{x, t; y, s) constructed in the previous 
section and estimate it, we need the following theorem. 

Theorem 5.1 (Lemmas 2 and 3 in [[T|). For each p > 0, let g{^', rj; p) be holomorphic function of 
(^', 77) in L"^^ c C""' X Cfor some /i > where 

L;-' = {(f , 77) 6 C""' xC;Imr]< nQRe t]\ + \Re ) - yu'^l/m }. 
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We assume the following estimate for g(^',rj;p). That is, there exist some constants C > and 
c > such that 

\g{e,mp)\ < Cd^'l + |77p^')'exp(-cp(|^'| + |77p/2))exp(C|Im^'|) (5.1) 
for {^', Tj) e L"-\ l<0,andp> 0. Then 

lv-'|2 1 ^2 

\Gix', t; p)\ < Cr^~5~^ exp [ - c ^ (5.2) 

where we set 

oo-iq 

G(x',t;p) = {2nr J e'^-^' J e'"^ gi^ , rj; p)d7ide 

R"-l -(xi-iq 

with an arbitrarily fixed positive number qfor p > and 

G(x',t;0) = lim G(x',t;p). 

Remark 5.2. 

(i) The theorem still holds even in the case the amplitude g depends on (jc„, j„). 

( ii) In /[T]/ the factor exp(C|Im^'|) does not exist in the estimate (|5.1I) . However, a slight modification 
of the proof given in [1] can include the case when the factor exists in the estimate. This was 
pointed out by Dr S. Nagayasu. 

We will apply Theorem 15.11 to estimate Y^^, in Appendix A. For this, we have to 
show that the assumptions of Theorem 15. II are satisfied in our case. To begin with we show that 
V^\xn,yn,^',r]), V^^{Xn,yn,^',r]), V\xn, y,,, ^' , are holomorphic in ,?/) e L^"' uniformly for 
(Xn,yn) with those Xn,yn satisfying the conditions attached to the definitions of V", V^^, V^. For 
simplicity we refer this property as uniform analyticity in L"^^ of V. Let y = (y,y) be either A or B. 
Then, the characteristic equation of each equation of (14.61) with respect to X has the form in terms 
of ^„ = iA 

Po^l + Pi(0^n + (P2(0-r) = 0, 

where po = -y„„, pi{^') = -2 ^"ij ynj^j, PiiO = - Yij^i^j and r = i?]. Then, the roots are 
given as ^„ = ~ ^ , where z~ = Jp^ - ^Poipi - "!") and +Im z~ > 0. By the ellipticity, for 

2po y 

some constant c' > 0, we have 

pI - ^PoPi < -c'\^'f (r e \ {0}, x: e U), 
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where U is an bounded open set in which PQ,pi,p2 are smooth. 

From the construction ofV^^, V^^, V^, the uniform analyticity of V in L^"' easily follows from 
the following lemma. 

Lemma 5.3. There exists yu > such that p\ - ApQP2 + 'Apoirj i [0, oo)for (^', rf) e L'^^^ and {Xn,yn) 
with those Xn,yn satisfying the condition attached to the definition ofV'. 



We will prove Lemma 153] by a contradiction argument. We first note that 

L;-i 3 (f , n)^Imn< K\Re nl + \Re ) - \Im (5.3) 

Suppose that for (^', 77) e L"^~^, there is m > such that p^ - 4pop2 + 4-poir] = m. Put a = -jnn, 
J3 = (j3u--- := -2(r„i, • • • ,r„„-i) and y = ini)i<i,j<n-i ■= i-7ij)i<i,j<n-i- Then, we have 

po = a<o, pi(r)=>e-r 

P2(0 = (rr)-r <0 for rer-i\{0} 
2 

m = 08 ® ;8) : ® - 4a(r^') ■ ^' + 4/a77. 
For simplicity, we denote = Re ^' and rjR = Re rj etc. Then, we have 

m = 08 ®;8) : (^;, ®^'j,)-(/3 ®P) : ® ^ + 2i(J3^/3) : (^;, 
- 4a{(r^;,) • ^'^ - (yO ■ - Sia(yO ' + ^iarj. 

That is, we have 

(5.4) 

08 <8>;8) : (^;, ® ^;) - 4a(y^'j,) ■ ^\ + la^R = 0. 
From the first equation, we have 

{-08 : i^'^ ® ^'^) + 4a(y^'^) ■ ^r} + {(fi®/3): ® ^ ' ^oc(yO ' ^'1) = - 4a/7/. (5-5) 

The left hand side (LHS) of (|5.5I) has the estimate LHS > c'|^^p-c"|^^p for some positive constants 
c' and c" . For the right hand side (RHS) of (15.51 ). by the definition of ^ we have from the second 
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equation in (15.41) 



RHS < -4ar]j < (-4a){/z(|7/«| + - m'' 
<nM{\^'^\' + \^f)-(-4a)fi-'\^f 



for some positive constant M. Thus, we have 



c(\eRf + < f^M(\e/ + (m > o). 



By taking /i > small, we have = = and hence t/r = by the second equation in (15.41) . 
Then, -Aarji = m gives rjj > 0. This contradicts to (15.31) . 

Next, we show the type of estimate (|5.1I) . Let V be one of V', / = 1 1, 12, 2 for either the case 
j„ > or yn < 0. Then, there exist positive constants c, C such that 



\V(Xn,yn,e,r])\ < C(iri + \r]\"^)-'cxp(-c\x,-y,m'\ + It?!^^')) exp (c|Imrik„ - (5.6) 



uniformly for iXn,yn) with those satisfying the condition attached to the definition of V'. 
Hence, by Theorem l5.1[ there exist a constant C > such that 



For the gradient estimate of r(x,t;y,s) we argue as follows. If we formally differentiate 
Y(x,t;y, s) by (1 < j < n - I), then the integrand is multiplied by i^j. Then, multiply what 
we got hy Xk - yki^ < k < n - I) and then integrate by parts with respect to ^i,. By these proce- 
dures, we end up with an integrand which satisfies the same type of estimate as (15.61) . Also, if we 
multiply by jc„ - y„ instead of multiplying by xj -yj(l < j < n - 1), we also have the same type 



of estimate as (15.61 ). because (jc„ - j„)(|^'| + \r]\^^^) exp ( - c\x„ - j„|(|^'| + \r]\^^^) is bounded. Hence, 



\Y{x, t; y,s)\< C{t - s) 2 exp ( - d 



,\x' -y'\^ + {xn-ynf 



) (t > s). 



t- s 



we have 



\V,,r(x, t; y, s)\ < C(t - s)-"^ exp ( - c 



\x-y\' 



) (t > s). 



t- s 



We can handle the derivative dx,X(x, t; y, s) in a similar way. Therefore, we have 



|V,r(.^, t;y, s)\ < C{t - sr"-i' exp ( - c'- —) (t > s). 



t- s ' 



(5.7) 
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